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INTRODUCTION 



Universal features of finite temperature and finite density variation of the QCD ground state realizing the SU{3)ii x 

I— I S'C/(3)l approximate chiral symmetry were understood with the help of the corresponding meson model p^. There were 

^S^ also numerous attempts to treat quantitatively the finite temperature restoration of chiral symmetry of the strongly 

^-1 interacting matter in the framework of this model. A common goal of these investigations is the description of the 

O , quark mass dependence of the nature of the finite temperature symmetry restoration :2 -4. . Another central issue is 

QJ the nature of the phase transition, when the baryonic density is varied. For its investigation one couples constituent 

T^ quarks carrying baryon number to the meson model [SHE]- 

Optimistically one could say, that the location of the characteristic points of the QCD phase diagram determined 
fvq with different variants of the model and in different approximations do agree with each other and with the results of the 
^ lattice field theoretical simulations within a factor of 2. In particular, improved agreement with lattice determinations 
CN of the QCD phase diagram were reported, when the Polyakov loop degree of freedom is coupled to the quark-meson 
OO model [5'-TI]. But the effective models are all strongly coupled, therefore one usually experiences large variations in 
^^ their predictions, when the simplest mean-field treatments are improved by taking into account quantum fluctuations 

. of the mesons and the constituent quarks |12| I13j. This circumstance limits the competitivity of their predictions. 

•/"i Some time ago we initiated the application of the expansion in the number of quark flavors for the description of the 

^^ chiral symmetry restoration in the two-flavor model 14, 15 . The approach was suggested by the SU{2) x SU{2) ^ 0(4) 

^—^ isomorphy and the fortunate fact that the N = oo solution of the 0{N) model is the textbook example of the application 

. . of the 1/N expansion [TS]. Recently, next-to- leading order (NLO) results were also presented for the pressure of 

J>" the relativistic 0{N) model and applied to the physical {N — 4) sigma-pion gas [T71 [T5|. In this context also the 

renormalizability of the NLO approximation was fully clarified and shown to be valid in various formulations of the 

J, model (with and without auxiliary field) ^IQji . 

d An analogous development for the three-flavor meson (and quark-meson) model is hindered, because even the leading 

order solution of the large-n limit of the U{n) x U{n) symmetric linear sigma-model is unknown. Progress would 

represent interest also for the Higgs sector of technicolor models of electroweak symmetry breaking [23 HI] ■ To our 

best knowledge, no published attempts exist, which would go beyond the partial large-n treatments of the 0(2n^)- 

symmetric nonlinearities. The results obtained with such an approach are questionable in light of the rather different 

finite temperature renormalization group behavior of the 0(2n^) and U{n) x U{n) symmetric models for n > 3 [Tll^. 

The goal of the present paper is to describe an approximate leading order (n = cx)) solution. Although it takes into 

account only two-loop contributions to the 2PI (two-particle irreducible) effective action, it definitely goes beyond the 

0(2n^) symmetric solution. It exploits in addition to the large n expansion an assumption concerning the mass spectra 

of the model. The region in the coupling space, where this assumption is valid can be estimated by determining the 

spectra from the approximate solution self-consistently. In this paper we present the construction of the renormalized 

version of this solution in some detail, and provide also an illustrative investigation of its range of validity. A deeper 

analysis of the mass spectra and the finite temperature features are left to a forthcoming publication. 
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The idea to impose an extra assumption on the mass spectra stems from the common practice of deahng with the 
4 scalars and 4 pseudoscalars defining the SU{2) x SU{2) symmetric meson model. One simply omits half of the 
fields, e.g. the 3 components of the scalar-isovector triplet and the pseudoscalar-isoscalar singlet with reference to their 
higher mass. One arrives this way at the Gell-Mann-Levy linear sigma model. We shall assume an analogous feature 
to occur in the spectra of the leading order solution at large n. We do not attempt the anomalous realization of the 
[/yi(l) symmetry, therefore in this limiting case all pseudoscalars will have the same mass. Based on this assumption, 
in a first approximation to the n = oo solution we retain only the quantum fluctuations of the light pion fields. 
After finding the propagators of the scalar fields, one calculates corrections to the pion propagator arising from the 
heavy scalar fluctuations. In principle, one can iterate this procedure until the solution of the large-n bare vertex 
approximation (BVA) is reached. It will be demonstrated, that the pion fields still obey Goldstone's theorem after the 
scalar corrections are included. We shall explore the divergence structure of the equation of state, the pion self-energy 
and the saddle-point equations (SPEs), and determine the counterterm pieces in the effective action necessary for their 
renor malizat ion . 

The paper is organized as follows. In Sec. II, the model is reformulated with two auxiliary fields. We elaborate also 
on the expected mixing structure of the solution. In Sec. Ill, the leading large-n Dyson-Schwinger (DS) equations are 
presented in BVA and the corresponding 2PI effective potential is constructed. In Sec. IV, we present an approximate 
solution, which is valid when the scalars are much heavier than the pseudoscalar fields of the model. Here also a 
simplified exploration of the range of validity of this assumption is given. The renormalizability of the proposed 
solution is investigated in Sec. V. In the divergence analysis we rely on a number of cutoff sensitive integrals collected 
in Appendix A. The paper ends with the conclusions (Sec. VI). 



II. FORMULATION OF THE MODEL WITH AUXILIARY FIELDS 



The Lagrangian of the model in its usual form reads as 
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L = Tr[9^M9^M^ - m^MM^] - ^(TrMM^)^ - ^Tr(MM^)2 + V2^/is°, 



2 .92. 



(1) 



where the algebra valued complex field M is parametrized with help of the generators of the U{n) group: 



M = (s°+i7r")r". 



ry rparph _A'^^ 



a ~ 0, 1, ..., n — 1. 



(2) 



In terms of the matrix elements, it can be written in the following form: 
where the U{n) vector t/° is defined as 



(3) 



U = —a \3 S -\- TT 7T ) — J STT, 



(4) 



with d"'"^ and /"''^ being the symmetric and antisymmetric structure constants of the U{n) group, respectively. Some 
useful relations involving them can be found in Appendix B., which are needed in order to obtain the subsequent 
equations. Two auxiliary fields are introduced by adding the following constraints to the Lagrangian: 



AL 



K-'^-/l^«'°>^ +<'•>'')' -K''°-/f^"^' 



In the sum C ^ L + AL the U{n) x U{n) — > U{n) symmetry breaking pattern corresponds to the shifts 

After the introduction of the auxiliary field variables and the shifts, the full Lagrangian has the following form: 
C 



(5) 



(6) 



1 r 
2 



{d^s^f + (df^TT^f - m'^hn'^v^ + 2V2^vs° + (s^f + {ir^f)] + V2^h{s° + V2n^v) 



V n 



(7) 



Next, we shortly describe the assumed structure of the solution and introduce some notations. The classical constraint 
equations show that the background v induces nonzero X, Y^ values and introduces mixing of the pair s*^, X and also of 
gO. ^Y°- for every value of the index a. We construct correspondingly a quantum solution, where the saddle-point values 
of X, y" are nonzero, F" = 0, a 7^ 0, and only the following 2-point functions do not vanish: 



[Gs°X, GsOyo, GxY°, Gs"s", Gyoyo, Gxx], 
[Gs'^s'',u=v:^Oj Gs^Y^,u=v^Oi Gy^Y^,u=v^o\j tr7r«'7r»',it=ii- 



(8) 



The sets in square brackets form mixing sets of fields: there is a 3-dimensional mixing sector, and there are n^ — 1 
identical copies of 2-dimensional mixing two-point functions. The n sector is diagonal. The equations presented in 
the next section show degeneracy of the 2-point functions for u ^ v ^ 0, therefore it is convenient to introduce the 
following short-hand notations: 



Ctt^tt" = G^, Gs^s"^ = Gs, Gs^Y'^ = GsY, Gy"Y^ = Gy, u ^ 0. 



(9) 



III. DYSON-SCHWINGER EQUATIONS AND THEIR TRUNCATION 
A. The equation of state and the saddle-point equations 



Standard rules 23 of constructing the derivatives of the quantum effective action F with respect to the fields leads 
to the following expressions: 



sr 

6X 
ST 

sr 



6T 



n 



(D + to2)5" _ ('mV2n2 _ 2iV^x)t;(5"° -K V2n^hd''° + i^^(Xs" + Gxs^) 



/Wi 



-(D -I- m^)^" + i^^^iXn'' + Gx^^) + i\ - d"''"(^''r'^ + G^.ya) ~ p''^{s''Y^ + G,.ya) 
ott" n y n I 



(10) 



Taking further functional derivatives of these expressions, one arrives at the equations of the 2-point functions (see next 
subsection). Also when one equates the derivatives of F to zero, the equations determining the physical value of the 
background v and the auxiliary fields X, y° are obtained (we do not distinguish the notations of the solutions of these 
equations from the variables taking arbitrary values). Using the Fourier representation of the propagators, from (10) 
one gets the two SPEs and the equation of state (EoS) (setting u = 0) at the leading order of the large n expansion as 



= —X + iny'2giV +i\ — n 



Gsik)+ / G,(fc) , 



-F° + iny/2i2V^ + iJ^n (J Gs{k) + J GAk)) , 



= -V2?i 



nv\ m —I- 



■ V^ 



lX-^^^-" 



Y") + V2nh + iy^n / Gsy(/c), 



(11) 



respectively. We note, that the equations for tt" and s" imply for u ^ vanishing expectation values automatically, 
taking into account the above assumed list of the nonvanishing 2-point functions. One should also note, that the 
auxiliary fields X, Y^ ^ 0{n). Comparing the solutions of the equations for the two auxiliary fields, one finds a relation 
valid both in renormalized and unrenormalized versions: 



V^y" = ^x. 

Also, it proves useful in the further discussion to introduce the combination 



(12) 



M^ = m^ - -(J2giX + ^232^ )• 
n 



(13) 



B. The propagator equations 



The second derivatives of the effective action yield the various 2-point functions. It is convenient to introduce the 
following tree-level inverse propagator commonly appearing in several two-point functions: 



iD^\p)^p'-M^ 



(14) 



In this section we substitute for all 3-point functions their classical expressions listed below in (151 and therefore close 
the coupled set of DS equations at the 2-point level (BVA). From (It]) one can read off the nonzero classical 3-point 
couplings: 



Y'^s^'s 



= r 



ya^t^c 



l92^abc!, 



^5''H{x-y)5{x-z), 



-d'""'S{x~-y)6{x~z), 



y.,.^. ^ -i^l^-^r'^Six - v)5{x - z). 
n 



(15) 



In the pseudoscalar sector the following DS-equations are found at leading order (the abbreviated notation / G[,-\G[,-\ 

/j, G'[.](fc)G[.](p + fc) is used): 



iG~^ - iD^' ~ ig2 J GviG, + G^) +ig2 J GsyGsy. 



(16) 



One notes the potential violation of Goldstone's theorem when comparing G„(fc — 0) or G-^o{k = 0) with the equation 
of state. Nevertheless, our approximate solution obeys this theorem due to a rather nontrivial relation between the 
relevant tadpole and bubble contributions (see next section). The previous relations also hint to a possible further 
dynamical violation of the remaining U{n) symmetry to SU{n) due to the coupling of the auxiliary variables Y"^ to the 
scalar fields s", which is proportional to the antisymmetric structure tensor f^^^ , This coupling can be seen explicitly 
in the structure of the V^ — s"" sector, where (if a ^ 0) one finds n? — \ identical mixing 2x2 DS equations: 



iG-^ - iD^^ - 192 J GsyGsy - i92 J Gy(Gs + G^), iGy' - -1 - ig2 J G^G^ - ^^ J{GsGs + G^G^), 



iG^Y = 2iy/g^v - ig2 / Gsy{Gs - G^). 
In the 3x3 mixing sector of {X, y°, s'^), the following 6 equations are obtained: 

iGx^ = -1 - igi / {GsGs + G-^Gtt), iG^ya = -iy/gig2 / {GgGs + G-^Gt,), 
*^7ox ^ 2i^grv - 2i^gig2 / GsGsY, iGyl = -1 - ig2 / {GsGs + G^^G^), 
iG~oyo = 2i^/g2V - 2i.g2 / GsGsY-, iG'^o = iD^^ - 2i.g2 / {GsGy + GsyGsy) 

One can collect the previously derived propagator equations and formulate them as matrices: 



(17) 



iQ 



{s,Y) 



iD^^ - Z.92 RGy{Gs + G,) + GsyGsy] 2i^v - tg2 J Gsy{Gs - G^) 
^i^v - tg2 J Gsy{Gs -G„) -l-if[A + J 2G,G,] 



(18) 



(19) 



-1 - igiA - iy/gig^A '2^^/glV - 2i^gig2 J GsGgy 

*^(x.y".sO) = ( -i^/9i92A -1 - ig2A 'i.i^v - 2ig2 J G^Gsy 

^2iy/g{v - 2iy/gig2 J GsGsY 2iy/g^v ~ 2ig2 J GsGsY iD^^ ~ 2ig2 J{GsGy + GsyGsy)^ 



Here we introduced the notation A = /(G^Gs + Gt^G-j^). 



(20) 



C. Construction of the two-loop 2PI effective potential 

In this subsection a 2PI effective potential {V^^^) is given, from which the previously determined equations can be 
obtained directly by functional differentiation. The functional depends on t h e ba ckgro und v, the auxiliary (composite) 
fields {X,Y^) and the 2-point functions. One should observe, that in (16 1, (191 and (20 1 the corrections to the tree- 






FIG. 1: Vacuum diagrams contributing to the leading order equation of G-^ (cf. the "setting sun" integrals of ( |23[ ) containing 
Gtt and proportional to v?). One pion line is cut corresponding to the variation of the 2PI effective potential with respect to the 
pion propagator. 



level propagators are given by 1-loop diagrams. Such terms arise from 2-loop vacuum diagrams contributing to the 
2PI effective potential by "cutting up" the lines corresponding to a propagator under functional differentiation (for 
an introduction see |24|). This means, that we have to draw 2-loop vacuum diagrams with the participation of the 
original and the auxiliary fields to reproduce the bubble contributions in the Dyson-Schwinger equations above (for the 
pion propagator equation the process is demonstrated on Fig. [I]) . This construction therefore corresponds to a specific 
truncation of the 2PI effective potential: from the infinite 2PI skeleton diagrams only a set of 2-loop graphs is taken 
into account. 

The following zeroth order (tree-level) propagator matrices from (l7| will also appear in the 2PI effective potential: 



iV 



1 _ / iDf^ ^ '2iy^v 



{'^,Y) 



2i^w 



- 1 







At this point it is convenient to introduce the rescaled variables 



-X, 



V 



f^y. 






(21) 



(22) 



which leads to M"^ = m? — i{x + y^). Using (11), (16), (19l and (20), one writes 



V 



2PI 



n^{m'^v^ - 2hv) + 



^ 2 , ^ l„a\2 \ ■2uP 1 ™\„,2 
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1 



G^{k)G^{k + p) + GMGs{k+p))GY{p) + 2GsY{k)GsY{k+p)Gsip) 
GY{k)Gs{k+p)-GsY{k)GsY{k+p))G^{p) 



.n 
+i-r92 



k J p 



n 
-«Tr52 



' k Jp 



1 

+^7;9i 



k J p 



Gs{k)Gs{k + p) + G^{k)GAk + p)]Gx{j>) 



+iV9mj J (Gsik)G,{k+p) + GAk)G^ik+p))GxYo{p) 
+i\92 I I {Gs{k)Gs{k+p) + G^{k)GAk+p))GYo{p) 



'S^VsIff^ / / GsY{k)Gs{k+p)Gs«x{p)+'^i92 / / G'.y(fc)G,(fc +p)G,oyo(p) 

J k Jp J k Jp 



+4^2 f J {G,ik)GY{k+p) + GsY{k)G,Y{k+p)]Gsoip) 



2 
-«ff2 



GYik)G^ik + p)G„o{p) + V, 



2PI 



k Jp 



(23) 



Since the tt and the s — Y sectors are of multiplicity (n? — 1), some 0(1) terms in ( [23| are not relevant for their leading 
order equations. These contributions were therefore not written out explicitly in ( |16p and (17). There are C(l) terms 
of the complete effective potential, which contribute only to the tt or s — F sector(at NLOjJ^ these are omitted from 



6 



(23 1. However, for the remaining propagators it is necessary to take into account the proper 0(1) terms, which are 



therefore included in (23 1. This means, that (23) is not a large-n expanded 2PI effective potential truncated at 2-foop 



level. Note, that the structure appearing in the fifth line corresponds to a setting-sun diagram with antisymmetrized 
vertex functions (e.g. ~ /"''^). 

In Vj^^ one collects all countcrterms, which should ensure the finiteness of the equations. One can introduce 
countcrterms to all independent pieces appearing in the mean-field part of the effective potential and independently 
of them also to the terms showing up in the inverse tree-level propagators. Note, that one might have to introduce 
counterterms also to pieces which are chosen to have a fixed numerical (zero or unity) renormalized coefficient. No 
counterterms are introduced corresponding to the setting-sun contributions, since these couplings are kept at their 
classical value. In this sense the most general form which we will need and is allowed by the structure of V"^^^ is the 
following: 



V^ 



2PI 



2 2 

4 ^ / ^\/9i92 



2^/5152 

J 

-(5yyG'y(fc) -f- 5xxGx{k) + (5yOyoGYo(fc) -|- 25xY°GxY°{k) 



[iSZ-j^k^ — Srri^ + iSxnX + iSy^y'^)GTr{k) + [SZgk^ — 5m1 + iS^sX + i5ysy'^)Gs{k)\ 



(24) 



(the self-energy corrections to G^y turn out to be finite). In the first line, countercouplings proportional to xy^,y'^,x 
will be also needed despite the fact that one chooses the corresponding renormalized values to vanish. Similarly, 
countercouplings in the second line belong to the contribution of the (pure) auxiliary propagators occurring in terms of 
the type TiV^^Q, in which all their renormalized coefficients are fixed to unity except the coefficient of Gxyo, which 
is zero. In the third line, countercouplings are introduced to the terms appearing in Dj^ and D~^. In the exact 
solution of the theory only unique quadratic and quartic counterterms should occur, but in any finite order of a 2PI 
approximation one has the freedom to choose the countercouplings in the counterterm functional independently pS] . 
For the determination of ( 24 ) one should analyze the divergence structure of the integrals appearing in the propagator 
equations. 



IV. A SELF-CONSISTENT ASSUMPTION ON MASS HIERARCHY 



The assumption, that the scalar sector is considerably heavier than the pionic simplifies the solution of the limiting 
form of the coupled Dyson- Schwinger equations valid at infinite n. In addition, we assume that the scalar masses are 
considerably larger than the amplitude of the symmetry breaking vacuum condensate. Each component of the propaga- 
tors in the mixed sectors will have a common denominator displaying the corresponding heavy mass. Therefore the only 
consequent way to neglect in a first approximation the heavy sector is to neglect the bubble contributions containing 
at least one component of Q(s,y) or Q(x,y'>^s2)- Then, in a first approximation only bubble diagrams exclusively built 
with G^ are included. This means that in (16) all bubbles are suppressed, while in (19) and (20) only the unique pure 
pion bubble remains. 



A. Retaining only the light pion dynamics 



The first consequence is that all n^ pion propagators are equal to their tree-level value, and therefore have the same 
mass: M^„ = Af^. The explicit form of the saddle-point equations can be written with the help of splitting the pion 
tadpole into finite (T^) and divergent parts [for the definitions of T^ ,T^ see Appendix A.] as 



J^ GAk) = T^ + rf + (Af2 - Ml)Tf\ 



and using the counterterms it becomes the following: 



= i[2v 



9i 
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t: 



1 



(^' '- (Af2 - Ml)Tf' + T/ - U.x + 5.y—=y 

' ^ V9i92 

i(2v' + rf ) + (Af2 - Afo^)Ti°' + Ti^) - (Syy^ + 5^y-^: 



2iSxQ 
-2i5. 



t/O- 



(25) 



(26) 



It turns out, that at this level the other counterterms (e.g. 5xv ,SxttiSxs, ^m ■, Sy^ ,Sys) are zero, theref ore for the sake 



of transparency we did not written them out explicitly in (26). Taking into account the definition (13) of M , the 



following counterterms renormalize both equations: 



<JX — Oy — Ij^ , 



5.0 -5y,^-\ (r(^) + (m^ ~ M2)rf )) , 5^y = ./gw2TT- (27) 



Because of the omission of the last ("heavy") term in the third equation of ( 11 ), the equation of state does not require 



any extra counterterm. The finite equations for the 1-point functions read as follows 
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The 2x2 propagator matrix of the {s,Y) sector simplifies to 



"^ - i{2v^ + T^), ^ - i{2v^ + T^), M\ = h. (28) 






(29) 



where we have introduced the function lT^(k^M) — ^i J Do{p)Dq{p + k) and also written down the contribution 
obtained from V^/"^ . By choosing 

-^rr^.frr, (30) 

one ensures the finiteness of the matrix elements. The squared scalar mass is determined by the zero of the determinant: 

1 

l-f/^(fc^Af(,,y),M)- 

The mass matrix of the {X, y°, s*') sector also becomes more transparent: 

1 + 2iSxx + giUik) 2iSxY0 + ^/9i92l7v{k) 2i^/g{v 



Ml.y) = M' + Ag,v' , o.rFn.^.. .,n ■ (31) 



^^(x,FO,s«) = 2«<5jfYO + y/gig2l-„{k) - 1 + 2i(5yoyo + ,g2/^(fc) 2i^/g2V ], (32) 

'^^\fg\'^ '^i^/giv iDq (fc) 



which after introducing the obvious counterterms 



6xx ^ ^giT^"\ 6xY'> = ^Vgi92T^°\ (^yoyo = -gaTJ"' (33) 

leads to a determinant equation completely analogous to the previous one: 

M^x,YO,sO) = M' + 4(gi + 92^- _^-^jF(1 m ^y (34) 

^ ' 1 - (.91 +52)J^(k = Ai(x,y",sO)7Ai') 

One can see, that the spectra resulting from the assumption we made for the mass hierarchy might be consistent with 
the outcome in the sense that the scalar sector (and the auxiliary fields hybridized with it) can be indeed heavier than 
the pionic one. 

In the next subsection we investigate the region of the parameter space, where the scalar masses become much heavier 
than the pseudoscalars ensuring the self-consistency of our approximate solution. 

B. Validity of the mass assumption 

First one has to note, that in the case when no explicit symmetry breaking term is added to the Lagrangian [h = 0) 
the mass assumption is true, since our approximation preserves Goldstone's theorem and therefore makes the pions 
massless, which means that they are "infinitely" lighter than the scalars. 

The interesting case is when h ^ 0. Let us define rs^y :— M(^gY)/M and rx,yo gO :— M(^x,y°.s°)/M . In order to 
obtain a proper region of the parameter space, we introduce a heaviness criterium: the scalars are heavy enough if 
relations Ts^y > fo^ '''x,y'>,sO > ^o hold simultaneously, where tq is a given number. It is somewhat arbitrary what value 
to choose for tq. In this explorato ry s tudy, we w ork with the convenient choice tq = 2, because the quantity I^ {k, M) 



appearing in both gap equations (ISlj) and (34) develops an imaginary part just for fc^ > 4M^ (two-pion threshold). 
Above the threshold the masses are defined as real parts of the complex solutions. Expressing v from the EoS in pS] ), 
the two relevant equations are the following: 



,2 ^ ,,2 , . ^ .-fSR4^fc = M(,,y),M) 



Kx) = M' + 432 ^ M ..^r^n.^..' :.M2 , r.:^ ;.,:_.. ^7^T2 ' (35) 







50 9 ^ 

FIG. 2: Region of the parameter space, where the mass assumption holds (using renormahzation scale Mo/^/jm^J = 2.5) 



M, 



(x,yo,.o) = A/f2 + 4(gi+52) 



1 - (.91 + g2m^ik = M(x.yo,,o), M) 



MMl-(.9i+.92)3?/^(fc = M, 



(x,yo,sO)i 



A^F + [(ffi+ff2)3/^(fc = M, 



(x.yo^gO) 



,M)P 



The bubble integrals are given by (A3). It is convenient to express all masses in proportion to the absolute value of the 
renormalized mass (-y/lm^l), which means practically to write the definition of the pion mass as M'^ = - 



IX ■ 



With the use of ( 28 1 and ( A4 ) , one gets 



M2 = -l + (.gi+.g2)( — 



IGtt- 



log 



eM2 



where Mq is the renormahzation scale. For fixed Mq, in the original units (36) determines M 



Wo- 



rn 



as a function of 



h/\m\'^,gi,g2- Plugging it into the scalar gap equations, they can be solved for M(^ y)/-^|r?T,2| and M(^x,y",s")/ \/\^\- 
Then, one can trace out the region, where the heaviness criterium is fulfilled. This region is the part of the positive 
9i — 52 octant (the stability region of the n — >■ oo theory) below the surface displayed in Fig. 2. As expected, the 
projection of the allowed region onto the gi—g2 plane shrinks for increasing value of h/\m\^. We have varied Mq/ ^/\m\^ 
in the interval (1, 10) and only a mild displacement of the allowed region could be observed without really changing the 
shape. This change can be balanced by an appropriate Renormahzation Group transformation of the quartic couplings, 
e.g. using g^ = g,(Mo). 

The modification of the allowed region occuring when the spectra is corrected by the fluctuations of the heavy scalars 
will be discussed in a forthcoming publication. 



C. Heavy scalar corrections 



Using the propagators of the coupled (s, F)-sector, one can start to systematically take into account the effect of the 
heavy degrees of freedom on the pion propagator, the EoS and the SPEs. In particular we now invoke the tadpole and 
bubble contributions to these equations evaluated with ( 29 ) . For this we write explicitly the components of the heavy 
scalar propagator matrix: 



G,ik) 



^(i-f4-(^)), GAk) = --^^ik^-An OMk)--^^2Vg-2v. 



d(fc) = (l-|/f(fc))(fc^-A/^)- 4.92^^2 



(37) 



With the help of these expressions, one readily writes the scalar tadpole contributions to the SPEs. The correction of 
the EoS is of imminent interest, since it has an important role in the discussion of the validity of Goldstone's theorem. 
Using the expression of GsY, one finds the unrenormalized equation of state: 



M^ + 2g2 



i h 

p d{p) V ' 



(38) 



Now we proceed with the unrenormalized pion propagators: 



iG-i(fc) = e-M^~2g2 



p2 - Af 2 



d{p) {k - p)2 - A-P ' 



J V 



p' - M' 



p d{p) {k - pY ~ NP 



d{p)d{k — p) Jp d{p)d{k — p) 

When one sets A; = 0, both propagator equations turn out to be the same. Comparing to the equation of state one finds, 
that the approximation fulfills the Goldstone theorem characterizing the U{n) x U{n) — >■ U{n) symmetry breaking. It 
will be demonstrated in the next section, that both equations receive the same counterterm contributions. 

Since the tadpoles of Gg and GsY play a role in the EoS and the SPEs, it is worthwhile to investigate their corrections, 
although these contributions should be considered only as NNLO "heavy" corrections. Substituting the leading order 
propagators into the second term on the right hand side of the equation of G~y (17 1, one finds for the E^y self-energy: 

^.rip) . .,. I GMp - k) (G.(fc) - GAk)) = -8,,V5..3 ^ ^^^^^^^^l^^^,^^^,^ . (40) 

which does not induce any counterterm, since it is finite. The heavy correction to G^^ is very similar to that of G^^ 
in (|39l): 

p2 - i\/2 



iG-^{k) = k^-M'^-g2 f 



p d{p) (fc - pY - NP 
^9lv^ I ZJm?, ^ - ^92 I ^" '" i:\.r^X' ""^' ^ (41) 



. r{p^~AP)il-g2l^{k-p)/2) 



Ip d{p)d{k - p) Jp d{p)d{k - p) 

The scalar corrections are taken into account also in the SPEs. They appear partly directly via the scalar tadpole. 



and also by the scalar correction of the pion tadpole. The structure of the two SPE's in ( 11 ) is identical in view of ( 12 ), 
therefore it is sufficient to investigate the equation, which determines x. The unrenormalized form of the saddle-point 
equation reads as follows: 

0^ — -2iv^-i[{Gs{k) + G.,{k)). (42) 

51 Jk 

The expression of the scalar tadpole is readily written down. The pion tadpole is expanded to linear order in the S^ 
self-energy contribution: 

-./^G.(fc).-.^I.o(^)+^^^,-^E., (43) 

where 

p2_M2 Ag,v^ , {p^-AP){l-g2l^{k-p)/2) 



E^(A:) =52 / 
Jp 



p d{p) 



{k-pY-AP d{k^p) d{k~p) 



(44) 



All these equations need (resummed) renormalization, which is discussed in detail in the next section. Note, that 
in this process the counterterms ( |27[ ), which were determined on the level of pure pionic fluctuations receive fur- 
ther contributions. Also, the yet unused counterterms (except SZg^SZ^^ counterterms proportional to wave function 
renormalization) will get nonzero values. 

V. DIVERGENCE ANALYSIS AND RENORMALIZABILITY 

The separation of the divergences in the EoS, the corrections to the pion propagator and the SPEs rely very heavily 
on our previous analysis of the NLO renormalization of the 0{N) model [15]. Most of the divergent integrals occurring 
in the present analysis can be made finite with the subtraction of the appropriate combinations of divergent integrals 
defined there. For the reader's convenience we list those, which arc used here in Appendix A. 

Let us start the divergence analysis with the EoS. The integral appearing in it can be expanded for large momenta 
in powers of 4(72 w'^: 

^■'^^"' ■ ^ (45) 



Ip dip) Jp (p2 - Af2)(l - g,I^{p)/2) V (p2 - M2)(l - g,I^{p)/2) 
Only the first term of this expansion is divergent and its divergence is given in (|A5[) with 7 = 52/2: 



p dip) 



= rP-^(M2-Mo2)Ty). (46) 
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From here one finds for the corresponding countercouplings of (24) 



Sml = -2g, (t(2) - ^(m^ - M',)T(')^ , S^^ = <S,„ = Sg^T^^^ 



(47) 



Next we proceed with the E^r s elf- energy. The second of the three bubble integrals appearing in the expression of the 
heavy scalar corrections to G^^ (391 is finite. The first integral can be written in power series with respect to 'ig2V^ as 



-W2 



1 



P 



AP 



d{p) {k - py2 - AP 
1 



'«32 



?2l' 



^ (p2 _ M2)(i _ g,lF^p)/2) V ' (p2 - Af2)(l - .g2/^(p)/2) 



AP 



{k - pY ~ A'P ■ 



(48) 



One recognizes again, that only the integral of the first term of the expansion is divergent. This is exactly the integral 
which was shown in Eq. (A7) of |19j to not have momentum dependent divergence. The same analysis leads also in 
case of the third integral to the same divergent piece. Therefore one concludes dZ^, = 0. 

In view of this, one can put fc = in the pion propagator and find its divergences. Since Goldstone's theorem is 
obeyed, the same divergences occur as in the EoS. As a consequence one has in (24) 



(5m^ = 5ml, ^x 



^yiT Oyy . 



(49) 



Using these results, one can find promptly the counterterms proportional to the G^ tadpole in (24). Since the divergent 
pieces of E^ coincide with those appearing in E^r, one has 



SZs = 0, 



Srn,. 



5ml, 



^XS 



Oyg Oyy . 



(50) 



These equalities reflect the fact, that the ultraviolet behavior of the tt and s multiplets are the same, irrespective of 
the symmetry breaking pattern and the structure of the mass spectrum. 



At this stage we have specified the expressions of the countercouplings of the counterterm functional ( 24 ) , up to 



the scalar corrections of the purely X, K'^-dependent divergences. The renormalized SPE of x takes a very transparent 
form, when the results of the above analysis are taken into account: 



0=^ 
91 



2« ( t-' + :^ /(G,s(fc) + G,(fc)) ) (1 + -5,.) + 5,x - 2i5^o + ^^ 
2 Jk ) V9i92 



SxyV 



(51) 



The last three terms represent the counterterms, which cancel those divergences of the tadpole integrals, which depend 
solely on x and y°. The consistency with the previous steps of the renormalization requires the cancellation of the 
"dangerous" product 5xv{'2v'^ + ^ Gs\f + S ^■^\f) — ^xvi^v"^ + T^ + T^) without the need for any new counterterm 
(which would react back on the EoS and the equation of the pion's 2-point function) . The unique common divergent 
coefRcient in front of this sum is expected to prevail in the exact renormalized equation, since it is the result of the 
Goldstone theorem and of the U{n) x U{n) symmetry valid in the ultraviolet regime. 

The contributions cancelling the above counterterm contribution come from the heavy corrections to the pion tadpole 
and from the scalar tadpole eventually also modified by the heavy corrections. This cancellation can not be complete 
on the level of the first round of heavy corrections, since then the scalar tadpole does not give any divergent term 
proportional to Tf which would be able to cancel SxvTf . However, this latter should be invoked first in the second 
iteration of the heavy corrections, therefore one simply omits it at this stage of the iteration, and no imbalance is 
detected in the divergences proportional to Tf . 

The actual form of the SPE (51 1 is the following after the first round of the heavy corrections: 



0^^~2i{v^ + ^ j{Gs{k) + GAk))] - 2i U^ + i /" Da{k)\ <5,, + ^x - 2i5xo + -^^Sxyv"" 



91 



V9192 



(52) 



Here one has to consider the pion tadpole as the sum displayed in ( 43 1 . The contribution of the scalar tadpole when 
expanded in powers of ^g2p can be written as follows: 



i / Gs{k) = -I 



. fc2 - Af 2 



1 



1 



72 w 



(fc2-M2)(l-.g24^(fc)/2) 



(53) 



Only these two terms, written out explicitly, contain divergences, therefore 



Gsik) 



div 



.[Ti^) 



<T^ 



(2) 



{aP 
{aP 



Mo')Tr] 






M2)2(l-g2/^(A:)/2) 



div 



(54) 
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Every term, except the last one on the right-hand side can be canceUed by appropriate extensions of the counterterms 
at most quadratic in x and 2/° [see (27l]. The heavy correction of the pion tadpole has the following expression: 



-i] G^{k) 



-W2 






?2W 



/p d{p) \{k- pY - K'P d{k - p) 
The part of the divergence which is independent of u^ is separated as 

1 



{p^~M^){l~g,I^{k-p)/2) \ 
d{k-p) ) 



(55) 



finite 



-ij G^{k) 



2ig2 



v—O.div 



fe (fc2 ^ M2)2 V7p {{k-pr - AP)il - g2l^{p,M)/2) 



finite_ 



(56) 



div 



Using (A12) of Appendix A., one finds 



A 



-2 I G^ik) 



= -2z52(Ad,„ + Brf„Af2) + 3iglT^iM)Ti'^ 



(57) 



v—Q, div 



The last term when substituted into (52) exactly compensates the dangerous term —iSxvT^{M), which arises from 
the —idxv J Dq term. The appearing new divergences (proportional to A^iy and B^iy) can be compensated by adding 
further heavy corrections to (27). The consistent cancellation of the divergences proportional to T^ could have been 
anticipated, since the counterterm proportional to the scalar tadpole left out from (51) at this level cannot lead to any 
divergence proportional to T^ . Therefore in this respect (52) does not differ from the exact relation (51). 



Next, one turns to the prospectively divergent terms proportional to v of the SPE contributed by the pion tadpole: 



, (fc2 _ Af2)2 

2 



; (p2 _ AP){{k-pr - M2)(l - g2l^{p)l2Y 



{{k-pf - A/2)(l - g2l^{k-p)l2){l - g2l^{p)/2) \ik-p)^ - AP p2 _ M^ 



(58) 



One finds by inspection, that there is no subdivergence in this expression. The first term of the square bracket produces 
an overall divergent piece 8ig2V^Ta . It can be found following the previous line of analysis, e.g. after changing the 
order of integration performing the /c-integration first. A careful, but rather lengthy analysis leads to the conclusion, 
that the second double integral is finite. 

Eventually one finds that in the sum the three divergent contributions proportional to w^ (Jq not annihilate. It 
is easy to see, at least partially, the source of this imbalance. When in a subsequent iteration round one uses in 
the scalar tadpole a propagator involving heavy corrections, it also produces a divergent contribution proportional to 
ig2V^Ta ■ This shows that in case of this single type of divergence for the cancellation one has to take into account 
contributions belonging to different iteration levels of the heavy contributions. Actually, all this is not unexpected, 
since the ultraviolet features of the integrals are independent of any hierarchy in the spectra. Therefore in practice, 
one renormalizes the saddle-point equations subtracting the remainder of the divergences proportional to u^ by hand. 



VI. CONCLUSIONS AND OUTLOOK 



A leading large-n solution was presented in a ground state, which breaks the U{n) x U{n) symmetry of the Lagrangian 
to U(n). In the construction of the solution, a light pseudoscalar/heavy scalar hierarchy of the spectra was assumed. 
The renormalizability of the equation of state and propagator equations, is a precondition for the investigation of the 
consistency of this additional assumption. This feature was demonstrated by constructing the counterterms to the 
equations above. The consistency of these counterterms required an additional subtraction for the renormalization of 
the saddle-point equations. It was shown, that the proposed solution explicitly fulfills Goldstone's theorem. The pions 
are generically light, therefore the assumed mass hierarchy will not be erased by heavy radiative corrections, as was 
demonstrated by our exploratory investigation. It is probably present in a large part of the (51,52) coupling-plane. 

The proposed procedure of constructing a leading large-n solution for the U{n) x U{n) — >■ U{n) symmetry breaking 
pattern can be developed further in several directions. In the range of validity of the mass hierarchy we shall study the 
finite temperature features of the proposed solution. The relevance of it would be largely strengthened, if a first order 
symmetry restoring transition were found in agreement with the renormalization group argument. For the applications 
to strong interaction phenomena at n = 3 one ought to introduce also the Ua{^) breaking effective (determinant) term 
as a sort of perturbation to this solution. A simple realization could be to include its contribution into the two-loop 
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2PI effective potential of the model and evaluate it with approximate large-n propagators constructed in this paper. 
Finally, we note, that one can make use of these propagators also in models, where constituent quarks are coupled to 
the mesons. 
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Appendix A. DIVERGENCES OF SOME RELEVANT INTEGRALS 

The divergences of the integrals listed below all can be read in somewhat scattered way in [19j . Here we summarize 
them for the reader's convenience. The divergences are expressed in terms of the following divergent integrals: 



t: 



(2) _ 



k^-M^ 



^d — 



(fc2 - Af2)2 



2^(2) 



Ak^-M^){l-^li{k)y 



rp(l 



, (fc2-M2)2(l-7/o^(fc))2 



l!{k), 



(Al) 



where Mq is an arbitrary normalization scale which makes these integrals infrared safe, 7 is a parameter which in case 
of our model equals 52/2, and I^(k, M) is the finite part of the bubble integral 



or in more explicit terms with real and imaginary parts separated: 

1 



-t: 



(0) 



m^ik.M) 



1 M2 



4M^ 
fc2 



-2./^^-larctan[l/Jl 



yi^^log 



I 



^A 



#+1' 



^,tF( 



1 



A]\,r2 



4M^1 
fe2 J' 



fc2 > 4M2 

4Af 2 > k^ >0 

0> fc2 



(A2) 



(A3) 



The finite part of the tadpole integral is 



T^{M) 



2\rpia) 



, P 



-J^-iM'-M^)Tr-T', 



.(2) 



AP . AP AP - M^ 



IQiT- 



'""^'AP 



167r2 



(A4) 



The divergences of the integrals below were found by replacing the mass parameter M sequentially (with help of 
subtractions and additions) by the normalization scale Mq: 



h = 



, (fc2_Af2)(l-7/f(A;,Af)) 



div 



= ri^) _ 3^(^22 _ Af2)ya) ^ f^„(M2), 



(A5) 



where I^{k, M) is the bubble integral defined with mass M. We need in our analysis also the slightly modified integral: 



ha^ 



1^ (fc2 - 4Af2)(l - 7/^(fc, M)) 
There is a logarithmically divergent integral: 



dii 



= fd,y{AP)-3-/ApTy\ 



k (A:2-M2)2(l-7/f^(fc,M)) 



T(I) 



div 



-7^^ 



(A6) 



(A7) 



The most challenging is the separation of the setting-sun integral, where one subdivergence is already explicitly 
found: 



fe (P - M2)2 yj^ {{k-pY - AP){l-^lF{p,M))J J 



(A8) 
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Using /i introduced above one can write 



, (P - M2)2 V7p ((fc - p)2 - M2)(l - 7/^(p, M)) 



+ ifd,y{M^) 



(A9) 



div 



After changing the order of the integrals, using (A19) of [19 , one is led gradually to the final form of its divergences 



h = - 



1/ ' 



d^v ^ U ^ ^*^' ^"^^ ^ 8^2 j y^ (fc2 „ 4M2)(1 - 7/i='(fc, M)) 



k (fc2 - M2)2 

T^^^ + (4M2 - Mo2)Ti°) 



^-^>'^^) + i 



(fd„ - 37M2Ti^)) 



-iT, 



k.(M2)(r]°)+/^(0,M) 



(AlO) 



The terms proportional to I^{0, M) can be rewritten with the help of a relation between the finite part of the tadpole 
integral T^{M) and the bubble integral at zero external momentum: 



M^l'"{0, M) = T^{M) + -7?-^(Af2 - Af2). 

Idtt" 



(All) 



In this way one finds a piece linear in M"^ with somewhat complicated looking divergent coefficients and a dangerous 
term resulting from an uncancelled subdivergence: 



/3 = i{Ad^. + BmvM^ - 3-fT^{M)Ty^), 



where 



A 



div 



T^^ - MlT^P - T, 



B, 



7 



div 


+ 


16^2 T^ a 


Ml- 


rp{0 
^d 


+ 


["^+16712 J 


^ a ■ 



Stt 



2^^d 



(0) 



^ divi 



(A12) 



(A13) 



Appendix B. USEFUL RELATIONS INVOLVING THE U{n) STRUCTURE CONSTANTS 

In the following relations, the summation over index e (where it appears) goes from to n2 — 1. 



foab — 0, doa 



b = \l —Oab: 

n 



n--l nr n -1 ,„ 

^ dan = (n^ - 1)W -, ^ dua = 5f)a(n^ - 1' 



i=l 



2 1 2 1 

n — 1 n — 1 

Y^ duj =0, ^ d, 

i=l,j/0 



aaf = y^n^^fO, 



'J 
/ ^ daaedecd = "^nScd, 2_^ du^dpcd = —{n ~ \)^cdi 



a=0 



a=0 



/ ^ ^ijc^ijd — ^ 



ri^-1 



1 + ScoSdQ ^(2 - (500(5^0) 



(^cd, 



/ ^ ^aic^aid ^ 



1 5" + ^cOf^dO 



(^cd, 



i=l a=0 



n^-1 



^ dacedade = n{l + ^c0^do)'^cd, ^ ficeftde = '^(1 - '^cO'^do)'^, 



cd- 



a=0 



(Bl) 
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